from scipy import log10, logspace, floor, ceil, pi, array
from numpy import concatenate

def freqint(sys):
    q = concatenate((sys.poles,sys.zeros))
    lowest = floor(log10(min(abs(q))))
    highest = ceil(log10(max(abs(q))))
    return logspace(lowest-2,highest+2,500)
